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A simple model is presented where the disappearance of domain 
walls and the associated production of "Fermi balls", which have been 
proposed as candidates for cold dark matter, are features which arise 
rather naturally in response to softly broken supersymmetry. 
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I. INTRODUCTION 

The possibility exists that the early universe experienced a sequence of symmetry 
breaking phase transitions, which could have resulted in the production of defects, 
such as monopoles, cosmic strings, or domain walls Furthermore, it is possible 



[ that supersymmetry could have been physically realized during early epochs, becom- 

ing broken after defect formation, at a lower energy scale. It is therefore quite natural 
to investigate models of defects within a supersymmetric context. Here, attention is 
focused upon a simple supersymmetric model constructed from a single chiral super- 
multiplet, which admits a domain wall solution interpolating between two distinct, 
but energetically degenerate, supersymmetric vacuum states. The domain wall for- 
mation arises from an exact discrete symmetry which is spontaneously broken. The 
initial supersymmetry of the model couples the fermion fields to the scalar fields in 
a prescribed way, and it is found that, as a result of this coupling, a fermion zero 
mode |4| forms within the core of the domain wall, where the fermion is essentially 
massless. Bound states can also exist which describe scalar bosons attached to the 
wall. 

The breaking of supersymmetry at lower energies can be described by the inclusion 
of soft supersymmetry breaking terms in the scalar potential. However, it is found 
that when the soft supersymmetry breaking terms are added to the Lagrangian, the 
exact discrete symmetry responsible for the formation of the stable domain wall is 
transformed into an approximate, or biased, discrete symmetry. This approximate 
discrete symmetry results in a domain wall network where each wall interpolates be- 
tween two different, energetically nondegenerate (and nonsupersymmetric) vacuum 
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states - a true vacuum state and a higher energy density false vacuum state [§,0. 
Thus, the terms that are added to the Lagrangian to break the supersymmetry also 
explicitly break the exact discrete symmetry. The model then resembles one recently 
proposed by Macpherson and Campbell (MC) wherein a biased discrete symme- 
try breaking results in the production of "Fermi balls" - tiny bags of false vacuum 
that are inhabited by a stabilizing Fermi gas. The Fermi balls emerge as an end 
product of the collapse and fragmentation of domain walls enclosing false vacuum 
protodomains. Massive, electrically neutral Fermi balls can be considered as candi- 
dates for cold dark matter. The simple model presented here therefore connects, in 
a rather natural way, the transformation of domain walls into Fermi balls and the 
breaking of supersymmetry. 

The supersymmetric model is presented in the next section, where a domain wall 
solution is found. Upon examining the response of the fermion and boson fields to 
the domain wall background, it is seen that a fermion zero mode forms inside the 
wall, and that there can exist bound states describing the attachment of bosons to 
the wall. The fermion mass vanishes in the core of the wall, and fermions near the 
wall can experience a short ranged, but strong, attractive force toward the core of the 
wall where it is energetically more favorable for them to reside. This effect will be 
of importance in the consideration of Fermi balls, since the wall will tend to absorb 
fermions to become populated by a Fermi gas of effectively massless fermions that 
can contribute a fermion degeneracy pressure. The soft supersymmetry breaking 
terms that are to be added to the Lagrangian are given in sec. III. It is easily 
seen that the inclusion of these terms causes the previous exact discrete symmetry 
associated with the formation of the domain wall to be exchanged for an approximate, 
biased discrete symmetry. The basic mechanisms proposed by MC relating to the 
disappearance of the domain walls through the formation of false vacuum bags, the 
collapse and fragmentation of the vacuum bags, and the production of Fermi balls 
are briefly reviewed. It is seen that the model presented here, with the breaking 
of supersymmetry, closely resembles the model presented for Fermi balls, allowing 
an inference that, at least in the context of the simple model presented here, the 
production of Fermi balls can arise from the breaking of supersymmetry and a discrete 
symmetry. A short summary forms sec. IV. 

II. THE SUPERSYMMETRIC MODEL 

Consider a supersymmetric model constructed from a single chiral superfield $ 
with component fields (cj),ip,F), where F represents the auxiliary boson field. The 
boson fields and F are complex scalar fields and the fermion field ip is a Weyl two- 
spinor. Let us write the scalar field in the form = A + iB, where A and B are 
real-valued. The superfield $ has a superspace representation |||| given by 

6) = <j>{y) + V29^{y) + 9 2 F(y), (1) 
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where = x» + iQuH and 6 2 = 66 = 6 a 6 a , a = 1,2. (We also have 0^ = 6 a ip aj 
6ip = 9a$ a , a = 1,2.) A metric g^ v with signature (+,—,—,—) is used. (See the 
Appendix for a brief description of the conventions and gamma matrices.) One can 
also define a Majorana 4-spinor \1/ in terms of the Weyl 2-spinors: 

a = 1,2, a = 1,2. (2) 



A. Lagrangian 



In terms of the chiral superfield the Lagrangian can be written as 

L = ($$)U- 2 + W($)\ e 2 + W{$)\ m , (3) 

where $ = $*, i.e., the "bar" and "star" symbols mean complex conjugation, Vl/($) 
is the superpotential, which will be defined shortly, and X\q2 stands for the 6 2 part of 
X, etc. By eliminating the auxiliary field F, the Lagrangian can be written in terms 
of the component fields as 

(4) 



L^ + L^ + Ly- V, 



where 



Lf = d^d^ = d^Ad^A + d^Bd^B, <j) = A + iB, 



T F _ i 
L K — 2 



'd 2 W\ , , (d 2 wY T T 



V 



dW 
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'dW\* 



(5) 
(6) 

(7) 



B. Superpotential and Scalar Potential 

To get a domain wall solution, let us choose the superpotential 

W = A$ (i$ 2 - a 2 ) (9) 

so that the auxiliary field is described by F* = —\(<f) 2 — a 2 ) , where a is a constant. 
From (|3|) the scalar potential is then given by 

V = F*F = A 2 (0 2 - a 2 )(4> 2 - a 2 ) = A 2 [(0</>) 2 - a 2 (0 2 + 2 ) + a 4 ]. (10) 

The scalar potential V = \F\ 2 > has minima located at F = 0, which implies 
that the (supersymmetric) vacuum states of the theory are located by <fi = ±a, i.e., 
the vacuum states of the theory are described by A = ±a, B = 0. It is seen that 
supersymmetry is respected in the vacuum, where V = 0. 
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C. Field Equations 



Recalling that (f) = A + iB, along with the expression for the Majorana spinor 
^ given by (0), the Lagrangian, expressed in terms of the real scalar fields A and B 
and the Majorana 4-spinor can be written in the form given by (|4]) with (using 

d 2 W/d(f) 2 = 2\cj)) 



L = d^Ad^A + d^Bd^B + \^^d^> + iX AW + B^/^ 
-A 2 [(A 2 - a 2 ) 2 + 2B\A 2 + a 2 ) + B 4 } , 



where (if)if) + ipip) = —i^ty and (ipip — ipip) = —^75^ have been used. With L = 
Lk + Lk + Ly — V, the field equations for A, B, and ^ follow from 



9* 9* 9* 1 ; 
Therefore, from ([TT| ) - ([14]) , the field equations for A, 5, and \I/ are given by 

□A + 2X 2 A(A 2 + 5 2 - a 2 ) - i\ M = 0, (15) 

UB + 2A 2 5(A 2 + B 2 + a 2 ) - i\ ^75^ = 0, (16) 



+ 2A(A + 575)^ = 0, (17) 



where □ = d^d^. 



D. The Domain Wall and Particle Masses 

Let us consider the real bosonic sector of the model where the fermionic fields 
vanish and the scalar field is real-valued, i.e. ^ = 0, B — 0. As boundary conditions 
for the A field we take A(x = ±00) = ±a. Then fllBD gives 

DA + 2X 2 A(A 2 - a 2 ) = 0, (18) 

which has as a static solution 
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X I 

Awix) = atanh — , w = —. (19) 
w Xa 

The solution given by (|1^) describes an ordinary domain wall of thickness w = 
We can notice that the domain wall can be associated with the spontaneous breaking 
of an exact discrete Z2 symmetry describing invariance of the real bosonic Lagrangian 
under A — > —A. Using this domain wall solution as a background solution, the 
response of the fields \l/ and B can be examined. It will be seen that there is a fermionic 
zero mode within the domain wall, and that there are domain wall- 1? particle bound 
states. 

In the vacuum states we have A = ±a, B = 0, \1/ = 0. Therefore, in vacuum 
the B particle mass is determined to be = 2Xa, and for the Majorana fermion, 
L^y^ = iXa^f^f = ^mptyty, which gives mp = 2Xa = m#, which is expected from 
supersymmetry in the vacuum. (For the vacuum state A = —a, the mass eigen- 
state Weyl spinors must be redefined by a phase rotation, and the Majorana spinor 
undergoes a 75 rotation, \& — > 75^.) 

In the core of the domain wall, A — > 0, and we find mp = 0, = y2Xa, so that 
the mass of each particle decreases within the domain wall. On this basis, we see 
that the particles are attracted toward the wall with a force F ~ —dm(x)/dx. The 
existence of scalar bound states and spinor zero modes is consistent with this picture. 
For the Majorana fermion we have m,p(x) = 2XA(x) (for x > 0, e.g.) and therefore, 



cosh — 



-2 



, which, 



by QT9D, the force of attraction can be estimated to be F ~ — -4 
for a thin wall, can be quite large in magnitude (but of short range, rapidly vanishing 
outside the wall's surface). 



E. Fermionic Zero Mode 

1. Static Zero Mode 

Upon setting A = A w (x), B = 0, the field equation for \l/ reduces to 

Yd^ + 2XA w ty = 0. (20) 

For the gamma matrices we have {7^,7^} = —2g^ u , 7 1 = i ( ^ ), (7 1 ) 2 = 1. 

Let us first look for a static solution of the form \1/ = ^(x). Multiplying (ED) by 7 1 
gives 

= -2XA wl ^. (21) 
Let us now write the Majorana 4-spinor \1/ in terms of 2-spinors r] and X- ^ = ( ^ j ■ 
We then have 7^ = z ^ oO ( " )=* ( ) " Therefore ' 
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rj 



-2i\A w 



o"i X 



0"! 



1 

1 



X J V ~ air l J 

The equations for r\ and x can be decoupled by writing 



si, 



7/ = ZCTiX- 



Then, by (0) and (E3 



A solution is given by 

?7 = rexp — 2A / A w (x')dx' 
L jo 



cosh — 
if. 



-2CXi77 



(22) 



(23) 



(24) 



(25) 



where r is an arbitrary constant Weyl 2-spinor. 

The Majorana condition — — 7 2 ^I / * = (where is the charge conjugate of 
\E0 i.e. 



icr 2 ?7* 



(26) 



must also be satisfied. Upon comparing ( |24"D and (p6[), we have o^ff = — aiT], or 
&1&2V* = so that with the help of <J\<Ji = ia 3 , we get 77* = 20377. We must 
therefore require that r* = ia 3 r. 



2. Traveling Waves 

Let us now regard $ to be a function of x, z, and t, i.e., ^/(x,z,t) 
rj(x,z,t) \ . . . r . -1-2 



—i(j\r\{x, z, t) 



where 77(2, z, t) = t(z, t) 



cosh - 



Then (£(]) implies that 



(j°d + 73 5 3 ) 



r(z,t) 

— zo"ir(z, t) 



cosh 



.1; 



77.; 



-2 



(27) 



which is solved by 



(d -a 3 d 3 )r(z,t) =0. 
This can be seen by multiplying ( |27| ) by 7 1 and using ry°~f 3 



(28) 
so that 



a 3 
^0 -<r 3/ 

reduces to the set of equations (<9 — <T 3 d 3 )r = 0, and (<9 + a 3 d 3 )aiT = 0, and 
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the second equation is automatically solved when the first equation is solved. Then 
(P§|) can be written explicitly as 



(<9 -<9 3 )0 \(n(z,t) 
(d + d 3 ) i I r 2 (z,t) 



0. (29) 



This is solved by 

T 1 (z,t)=T 1 (z + t), T 2 (z,t)=T 2 (z-t). (30) 

Therefore, r can be written as 

*.fl-(?" +t) ) + (v«))- ^ 

so that r, and hence can contain a linear combination of "up" and "down" moving 
waves. 



F. B Particle Bound States 

Now let us set \I/ equal to zero and examine the B field in the domain wall 
background. From the field equation for B, we have 

UB + 2\ 2 B(A 2 W + B 2 + a 2 ) = 0. (32) 

Now linearize, and look at small fluctuations about B = to obtain 



UB + 2AV 1 + tanh^ — ) B — 0. (33) 
V wj 

Writing B(x, z, t) = b(x) sm(kz — ut + 5), (j33|) reduces to 

-d 2 x b + 2X 2 a 2 [ta,nh 2 (x/w)}b = E\ E 2 = oo 2 - (k 2 + 2AV). (34) 

This is a Schro dinger- like equation with an attractive potential that can accommodate 
one or more bound states |10| with < E < \^2\a. We therefore infer that real scalar 
B particles can be localized within or near the core of the domain wall. For E > \^2\a 
there can exist a set of states describing the scattering of B particles from the domain 
wall. 
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III. SOFT SUPERSYMMETRY BREAKING AND FERMI BALLS 



The supersymmetry that exists in the vacuum states of the above model can be 
broken by adding soft supersymmetry breaking terms to the scalar potential [fTI]. The 



types of soft terms allowed here include a scalar mass term of the form /x 2 (/></> and a 
trilinear scalar interaction of the form [W($)|e=o + c.c] = go[<f) 3 + <ft 3 }. Let us also 
add a (dynamically irrelevant) constant Vq and therefore define the potential term 

V B = /J ! 00 + g (0 3 + 3 ) + V . . 

= ^(A 2 + B 2 )+g (A 3 -3AB 2 )+V , 1 ' 

where the constant V can be used to adjust the vacuum energy of the true vacuum 
state to zero. The total potential can be written as 

Vi = V + V B . (36) 

Notice that not only has the original supersymmetry been broken, but the discrete 
Z2 symmetry associated with the reflection A — > —A has also been broken by Ve- 
in fact, the model now resembles the kind of model that was introduced for the 
description of a Fermi ball. Therefore, we have the possibility that a simple super- 
symmetric model, which possesses an exact discrete symmetry before the breaking of 
supersymmetry, can be left with an approximate discrete symmetry after the breaking 
of the supersymmetry, so that energetically nondegenerate vacuum states develop - 
a true vacuum state and a higher energy false vacuum state. Two different domains 
are separated by a domain wall, which can rapidly fragment into a mist of Fermi balls 
by the mechanism described by MacPherson and Campbell (MC). Massive neutral 
Fermi balls which couple only weakly with ordinary matter are then candidates for 
cold dark matter. Let us briefly review the Fermi ball model proposed by MC and 
then examine the model presented here to see how it can describe Fermi balls. 



A. Fermi Balls 

The basic scenario described by MC for the production of Fermi balls can be briefly 
summarized in the following. (For discussions of biased discrete symmetry breaking 
and domain walls, see also refs. f5|||.) Consider a simple model of a self-interacting 
scalar field ip and a Dirac fermion field ip which is strongly coupled to the scalar field 
(p. The system can be described as follows: first consider the Lagrangian 

L = U^d^-^-^lf. (37) 

L possesses a discrete Z 2 reflection symmetry, and the degenerate vacuum states 
of the theory are described by (p — ±y?o- A domain wall solution of the form ip = 
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<^ tanh^|J interpolates between the two distinct vacuum states, with tp — > ±y? as 
x — > ±00, and in the core of the domain wall </? — > 0. The surface tension is equal to 
the surface energy density, given by 

o = HM. (38) 

Let us now consider changing the exact discrete Z 2 symmetry to an approximate 
symmetry. The symmetry breaking results in the formation of two different, nonde- 
generate, vacuum states that form protodomains of true and false vacuum. Let the 
difference in the energy densities of the two vacuum states be represented by A. In 
this biased discrete symmetry breaking a domain wall can form which interpolates 
between the true and false vacuum protodomains. The exact Z% symmetry can be 
exchanged for an approximate symmetry by adding, for example, a Z2 symmetry 
breaking term A(<p) to the Lagrangian. The asymmetry that is introduced can re- 
sult in the formation of finite sized "false vacuum bags" - regions of false vacuum 
protodomain enclosed by domain wall |5|||. These vacuum bags can collapse, and 
result in the conversion of false vacuum into true vacuum and the disappearance of 
the domain walls. 

Let us now consider a Dirac fermion if> that is strongly coupled to the scalar field 
ip through a standard Yukawa coupling. We replace the Lagrangian Lq with 

L = $(yfdp + iGp)^ + \d^d^ - ^ (p 2 - tpl) 2 + A{<p). (39) 

[A relative factor of (— i) appears in the Yukawa term in (|39|) due to our choice of 
representation for the gamma matrices.] The fermion acquires different masses in the 
different protodomains due to its coupling with <p, but the fermion becomes effectively 
massless in the core of the domain wall where <p — > 0. It is therefore energetically 
favorable for the fermion to reside inside the wall, and it is assumed that fermions 
near the wall will become absorbed by the wall, so that the wall is quickly populated 
by massless fermions. As an estimate, we can think of a force of attraction acting on 
the fermions by the wall to be given roughly by f(x) ~ —dM(x)/dx ~ —Gd(f(x)/dx. 
In the thin wall approximation, we can then regard the domain wall as possessing a 
two dimensional Fermi gas of massless fermions. The Fermi gas pressure can have 
a stabilizing influence by counteracting the wall surface tension and false vacuum 
pressure contributions. For a spherical vacuum bag, the collapse can be halted when 
the bag has a radius R, which minimizes the total energy E, and is related to the 
number of fermions A" inhabiting the wall. 

However, a vacuum bag of energy E and radius R is not stable against flattening 
into a "pancake" shape. The tendency to flatten thus results in the fragmentation of 
the vacuum bag into many smaller ones. The fragmentation process halts when the 
thin wall approximation is no longer valid, i.e. when the typical radius of curvature 
of a bag becomes comparable to the wall thickness. In this limit, the configuration 
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is better represented by a "Fermi ball" which can be thought of as a tiny vacuum 
bag with essentially no false vacuum in the interior - a nontopological scalar field 
configuration consisting mostly of the domain wall inhabited by a three dimensional 
Fermi gas. By equating the minimum size of the stabilized Fermi ball R m m to the wall 
thickness 5, and assuming the stable Fermi ball to be spherical, the typical stabilizing 
radius at which the collapse and fragmentation process stops is estimated to be 

Rmin ~ T~ ' (40) 

Aipo 

The presence of a Fermi gas inside the domain wall is crucial to the formation of 
Fermi balls in this model. In order that fermions and antifermions do not undergo 
annihilation processes that leave no Fermi gas inside the wall, it is sufficient to assume 
that there is a net fermion antifermion asymmetry, so that annihilation processes 
which may occur will eventually stop when all of the antifermions (or fermions) have 
been consumed, leaving a Fermi gas of fermions (or antifermions). 

MC estimate that a Fermi ball would contain about 50 fermions and have a mass 
of roughly 100(^o GeV, where ipo is expressed in units of GeV. If the Fermi balls are 
constructed from a new Dirac fermion that has no standard model gauge charges, then 
the Fermi ball would likely be a neutral, heavy, nonrelativistic particle interacting only 
very weakly with ordinary matter. In this case, Fermi balls would form a candidate 
for cold dark matter. 



B. Supersymmetry Breaking and Fermi Balls 

In the model presented here, the soft supersymmetry breaking terms are embedded 
in the potential term Vb, given by fl35|). Thus, by ( [TT]) and ( |55D the total Lagrangian, 



written in terms of the real scalar fields A and B and the Majorana field is given 
by 

Lx = d^Ad^A + d^Bd^B + i^fytf + iA[A*tf + B^^l] 

-X 2 [(A 2 -a 2 ) 2 + 2B 2 (A 2 + a 2 )+B i ] ' (41) 

-fi 2 (A 2 + B 2 ) - g (A* - 3AB 2 ) - V . 



This model closely resembles that presented by MC, except that (f41|) contains an 
additional scalar field B and the fermion here is Majorana, rather than Dirac. [Let 
us also assume a parameter range that keeps the vacuum expectation value B vac = 

2 

after supersymmetry breaking. For example, we could require that ^-j << 



a 



2 



2# << a, and jj 2 — 3g a > 0, so that after the breaking of the supersymmetry 
Bvac — and the vacuum values of A are only slightly shifted from ±a.] Upon 
setting B = 0, (flip is seen to have the same form as (0). The supersymmetry 
breaking terms embedded in Vb also break the exact discrete symmetry associated 
with A — > —A in the supersymmetric version. We therefore expect the Fermi ball 
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scenario to be realized in this broken supersymmetric model, as well. Following 
the reasoning of MC, we expect a Fermi gas to remain within the domain wall after 
possible fermion antifermion annihilations if there is a fermion antifermion asymmetry, 
or when it becomes energetically unfavorable for massive particles, like the A and B 
scalar bosons, to be produced outside the wall. The Fermi balls associated with this 
model are neutral, and again can be considered as candidates for cold dark matter. It 
can also be pointed out that a slightly more complicated model |L2|], composed of two 
interacting chiral superfields (and hence two Majorana fermion fields, or equivalently, 
a Dirac field), that has the same basic features presented here might be implemented. 
The Majorana fermion can then be replaced with a Dirac fermion. It is conceivable 
that a model of this type could be constructed where the scalar and spinor fields 
are allowed to have standard model gauge couplings. But one of the points to be 
illuminated in this work is that the MC type of model, which leads to the prediction 
of Fermi balls, can emerge rather naturally in an initially supersymmetric theory 
where the supersymmetry gets broken, along with the exact discrete symmetry. The 
breaking of the exact discrete symmetry can remove the potentially hazardous domain 
wall problem, and also give rise to the production of Fermi balls. 



IV. SUMMARY 

Because there exists a strong possibility that (1) the early universe underwent 
a set of symmetry breaking phase transitions, during which defects may have been 
formed, and (2) supersymmetry was physically realized at the time of defect produc- 
tion and was broken at a somewhat later time, it becomes relevant to consider field 
theoretic models of defects within the context of supersymmetry. Here, a simple do- 
main wall model, constructed from a single chiral superfield, has been examined. It is 
found that the fermions become massless inside the domain wall, where a zero mode 
forms, and that there are bound states describing real scalar bosons attached to the 
wall. (These types of results have been examined previously for a supersymmetric 
model constructed from two interacting chiral superfields [T2J, but the effects of soft 



supersymmetry breaking terms are more easily examined in the single field model 
presented here.) 

An exact Z2 discrete symmetry in the bosonic sector of the theory, associated 
with the reflection symmetry A — > —A, gets broken explicitly to an approximate 
biased discrete symmetry when soft supersymmetry breaking terms are added to the 
Lagrangian. The model then closely resembles one describing "Fermi balls" , which are 
scalar field configurations stabilized by a Fermi gas exerting a degeneracy pressure. In 
the context of the simple model presented here, we therefore expect the production of 
domain walls, followed by a process wherein the walls form false vacuum bags (due to 
the transformation of the exact discrete symmetry into an approximate one), which 
collapse and fragment, finally resulting in the production of Fermi balls. Thus, in 
the model presented here, the production of Fermi balls is closely associated with the 
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breaking of supersymmetry. 



APPENDIX A: CONVENTIONS 

Some of the notations and conventions are briefly listed here. A metric g^ v is 
used with signature (+, — , — , — ). Aside from the metric, the notation, conventions, 
and gamma matrices used conform to those of ref. ||. The gamma matrices can be 
written in the form 

with 

<7* = (1,*), a» = {l,-a), (A2) 
where a represents the Pauli matrices. Then 



and 75 is given by 



75 = tVtV = M J _°i J ■ (A4) 



The gamma matrices have the properties 

{Y,Y} = -W, {7 M ,75} = 0, 75 = -75, (75) 2 = -l- (A5) 
A Majorana 4-spinor \1/ is expressed in terms of the Weyl 2-spinors if) and if> by ^ = 
( ^& j an d we use the summation conventions for Weyl spinors [with ip a = (if> a )*] 



0--C'- n . 0--0,c". a = 1,2, a = 1,2, (A6) 
with e metric tensors (for raising and lowering Weyl spinor indices) 

= (s^) = ia 2 , {e a p) = (e. ) = -ia 2 , e 12 = 1 = e i2 . (A7) 



[1] A. Vilenkin, Phys. Rep. 121, 263 (1985) 



12 



[2] A. Vilenkin and E.P.S. Shellard, Cosmic Strings and other Topological Defects (Cam- 
bridge University Press, 1994) 
[3] E.W. Kolb and M.S. Turner, The Early Universe (Addison- Wesley, 1990) 
[4] R. Jackiw and C. Rebbi, Phys. Rev. D 13, 3398 (1976) 
[5] G.B. Gelmini, M. Gleiser, and E.W. Kolb, Phys. Rev. D 39, 1558 (1989) 
[6] J.A. Frieman, G.B. Gelmini, M. Gleiser, and E. W. Kolb, Phys. Rev. Lett. 60, 2101 
(1988) 



[7] A.L. MacPherson and B.A. Campbell, Phys. Lett. B 347, 205-210 (1995); [Also, [hep- 



ph/9408387 Biased Discrete Symmetry Breaking and Fermi Balls. 



[8] See, for example, P. Srivastava, Super symmetry, Superfields and Supergravity: an In- 
troduction (Adam Hilger, 1986) and references therein. 

[9] See, for example, J. Wess and J. Bagger, Supersymmetry and Supergravity, Second 
Edition (Princeton University Press, 1992) 
[10] P.M. Morse and H. Feshbach, Methods of Theoretical Physics (McGraw-Hill, 1953) p. 
1651 

[11] See, for example, R.N. Mohapatra, Unification and Supersymmetry, 2nd ed. (Springer, 

Berlin, 1992) 
[12] J.R. Morris, Phys. Rev. D 52, 1096 (1995) 



13 



